Nonlinear localized waves in a periodic medium 
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We analyze the existence and stability of nonlinear localized waves in a periodic medium described 
by the Kronig-Penney model with a nonlinear defect. We demonstrate the existence of a novel type 
of stable nonlinear band-gap localized states, and also reveal an important physical mechanism of 
the oscillatory wave instabilities associated with the band-gap resonances. 
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Wave propagation in periodic media is associated with 
many interesting physical phenomena plB. Modern 
technology allows to create different kinds of macro- 
and mesoscopic periodic and layered structures such as 
semiconductor super lattices and heterostructures, mag- 
netic multilayers possessing the giant magnetoresistance, 
multiple-quantum- well structures, optical waveguide ar- 
rays, photonic band-gap materials and photonic crystal 
fibers, etc. The main feature of different periodic struc- 
tures (which follows from the classical Floquet-Bloch the- 
ory) is the existence of forbidden frequency band gaps (or 
stop bands) where linear waves undergo Bragg reflection 
from the periodic structure pL However, many of the re- 
cently fabricated periodic structures exhibit pronounced 
nonlinear properties that give rise to qualitatively new 
physical effects such as multistability of a finite nonlinear 
periodic medium || and energy localization in the form 
of gap solitons Q|. Such effects are usually analyzed in 
the framework of the coupled- mode theory [0 , and they 
are associated with the nonlinearity-dependent tuning of 
the stop band as the wave intensity is increased. 

Another fundamental reason for the wave localization 
is the presence of defects. As a matter of fact, artificially 
introduced inhomogeneities are often used as a power- 
ful means of controlling wave scattering. However, the 
defect-induced localization in nonlinear and periodic me- 
dia is largely an open area of research. We stress that 
the defect response is determined by the local field am- 
plitude, and therefore the standard averaging procedure 
can not be directly applied to this kind of problems. 

Recent experiments, e.g. the observation of optical 
gap solitons ^ and the control of coherent matter waves 
in optical lattices 0, as well as theoretical results such 
as the discovery of the oscillatory instability of gap soli- 
tons ||, call for a systematic analysis of nonlinear ef- 
fects in periodic structures and band-gap localized states 
beyond the approximation provided by the coupled-mode 
theory. Such an analysis is crucially important for deter- 
mining stability of nonlinear waves in periodic structures 
because the wave instabilities can appear due to the mode 
coupling to other bands. 

In this Letter, we present the first analysis of the exis- 
tence and stability of nonlinear localized waves in a peri- 
odic medium with multiple gaps in the transmission spec- 
trum that is valid beyond the coupled-mode theory. We 
consider a simple model where waves are localized in a 



layered medium by an intensity-dependent defect (or, in 
other words, they are guided by a thin-layer nonlinear 
waveguide). Assuming the applicability of our results to 
a variety of different physical systems (see below), we 
describe four qualitatively different situations and char- 
acterize, in the framework of a unified and systematic 
approach, the properties of two types of nonlinear local- 
ized waves P]. For the first time to our knowledge, we 
analyze stability of nonlinear localized waves in a periodic 
medium and reveal an important physical mechanism of 
wave instability associated with the band-gap resonances. 
We demonstrate also that several types of stable band- gap 
localized states can exist in the presence of nonlinearity. 

Model. We describe localized waves in superlattices 
using the nonlinear Schrodinger equation for the wave 
envelope 4>{. x i *)> 
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where I = \ip\ 2 is the wave intensity, t is time (or prop- 
agation variable), x is the spatial coordinate, and the 
real function x) describes both nonlinear and peri- 
odic properties of the medium. We note that the sys- 
tem (Q) is Hamiltonian, and for localized solutions the 
power P — f_™ I(x) dx is conserved. 

We seek stationary localized solutions of Eq. (|l]) in the 
form ip(x,t) = u(x)e iult , where ui is the normalized fre- 
quency (or the propagation constant, in optics), and the 
real function u{x) satisfies the equation: 
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We assume that the superlattice is linear, and nonlin- 
earity appears only through the properties of an em- 
bedded localized defect. Then, if the corresponding 
width of the wave envelope is much larger than that 
of the defect, the inhomogeneity can be modeled by a 
delta-function and, in the simplest case, we can write 
T{l\x) — v(x) +S(x)G(I), where the function G(I) char- 
acterises the properties of the defect, and v{x) = v(x+h) 
describes an effective potential of the superlattice with 
the spatial period h. 

For such a local nonlinearity, the localized waves can 
be constructed with the help of a matching condition, by 
using the solutions of Eq. with T(I;x) — v(x), pre- 
sented in the form of the Bloch-type wave functions pj . 
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Band-gap structure and localized waves. If the effective 
periodic potential v(x) is approximated by a piecewise- 
constant function (the so-called Kronig-Penney model), 
the solution on both sides of the defect can be decom- 
posed into a pair of counter-propagating waves with the 
amplitudes a(x) and b(x), 

u b {x) = a{x)e-^ x + b(x)e^ x)x , 

where pk{x) — yoT— v(x) is the local wavenumber. As 
follows from the Floquet-Bloch theory, for a Bloch-wave 
solution the reflection coefficient r(x) = b(x)/a(x) is a 
periodic function, i.e. r(x) = r(x + h), and it is found 
from the following eigenvalue problem: 

where T(x) is a transfer matrix that describes a change of 
the wave amplitudes {a, b} after one period (x, x+h). An 
explicit expression for the transfer matrix can be found 
in Ref. , where it was also proven that det T = 1 , and 
two linearly independent solutions of Eq. (^|) correspond 
to a pair of the eigenvalues r and r _1 . Relation t(uj) 
determines a band-gap structure of the superlattice spec- 
trum: the waves are propagating, if |r| = 1, and they are 
localized, if |t| ^ 1. In the latter case, a nonlinear defect 
can support nonlinear localized waves, and the wave am- 
plitude at the defect is determined from the continuity 
condition (lO) at x = [I = 1(0)]: 

Go = G(I ) = CM- (4) 

Here C = (C + + r)*=o, C* = M ± (l - r ± )(l + r*)' 1 , and 
we denote with "+" and "— " the lattice characteristics on 
the right and the left side of the nonlinear defect, respec- 
tively, i.e. v(x) = u + (\x\), for x > 0, and v(x) = v~ (\x\), 
for x < 0. Note that the reflection coefficients should 
correspond to the Bloch-wave solutions with the asymp- 
totics \u(x — > ±oo)| — > 0. 

Relation (||) allows us to identify different types of 
nonlinear localized states. We notice that such local- 
ized states, supported by an attractive nonlinear defect 
(Go > 0), exist in the so-called waveguiding regime, when 
C(w) > 0. Additionally, localization can occur at a repul- 
sive nonlinear defect (Go < 0) in the anti-waveguiding 
regime, provided £(^-0 < 0. A characteristic band-gap 
structure for a two-component superlattice is presented 
in Fig. |l| (top). Using the guided- wave terminology, we 
notice that the first (semi-infinite) band-gap corresponds 
to the conditions of the total internal reflection (IR) and, 
therefore, its dispersion properties should be similar to 
those of the conventional waveguiding regime. Localized 
waves in this band resemble conventional solitary waves 
modulated by a periodic structure [see Fig. [j](b)]. In 
contrast, both waveguiding and anti-waveguiding regimes 
can occur at smaller oj, where band-gaps appear due to 



the resonant Bragg-type reflection (BR) from the peri- 
odic structure, so that the localized waves are somewhat 
similar to the gap solitons composed of mutually coupled 
backward and forward propagating waves [see Fig. |](a)]. 
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FIG. 1. Top: power vs. frequency dependences for 
the localized waves: solid — stable, dashed — unstable, and 
dotted — oscillatory unstable. Middle: real (dotted) and 
imaginary (solid) parts of the eigenvalues associated with the 
wave instability. Shading marks "waveguiding" (white) and 
"anti-waveguiding" (dotted) localization regimes inside the 
band gaps. Bottom: the localized states corresponding to the 
marked points (a,b) in the top plot; shading marks the areas 
with smaller v. The lattice parameters are h — 1, v(x) = 
for n — 1/2 < x/h < n, and v(x) = 30 for n < x/h < n + 1/2, 
where n is integer. The impurity characteristics are a = 0.5, 
= 1. 

To study the linear stability, we consider the evolution 
of small-amplitude perturbations of the localized state 
presenting the solution in the form 

iP(x, t) = |n(x) + v{x)e^ t + w*(x)e^' t ^ e tut , 

and obtain the linear eigenvalue problem for v(x) 
and w(x). Then, from the condition of mode lo- 
calization we define the so-called Evans function, 

Y(i) = [G x -au + l)][Gx-au-l)\ - (Gi-Go) 2 , 
where Gi = Gq + IoG'(Io), the intensity Jo is calculated 
for an unperturbed solution, and the prime stands for the 
derivative. Eigenvalues are found as zeros of the Evans 
function, ^(7) = 0, and the corresponding eigenmode 
solutions fall into one of the following categories: (i) in- 
ternal modes with real eigenvalues describe periodic oscil- 
lations ("breathing") of the localized state, (ii) instability 
modes correspond to purely imaginary eigenvalues, and 
(iii) oscillatory instabilities can occur when the eigenval- 
ues are complex. 
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FIG. 2. An example of a resonance that occurs between 
an internal mode of the localized wave and a bang-gap edge 
and leads to an oscillatory instability. 

To demonstrate the basic stability results, we consider 
a localized defect possessing a cubic nonlinear response, 
G(I) — a + (31. Under proper scaling the absolute value 
of the nonlinear coefficient (3 can be normalized to unity, 
so that (3 — +1 corresponds to self-focusing and (3 = — 1 
to self-defocusing nonlinearity. Localization depends also 
on the sign of the linear coefficient a, which defines the 
defect response at small intensities: attractive if a > 0, 
and repulsive, otherwise. Therefore, below we consider 
four qualitatively different examples that correspond to 
the different signs of a and (3. 
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FIG. 3. Power vs. frequency, and (a-c) wave profiles for 
a = —5, (3 = 1. Notations are the same as in Fig. 

Self-focusing nonlinearity. First we consider the prop- 
erties of the localized waves supported by a defect with 
a self- focusing nonlinearity ((3 = +1), and attractive lin- 
ear response (a > 0). Such waves already exist in the 
linear limit at the frequencies u>b defined by the equal- 
ity Ci^b) — ct, and they correspond to the waveguid- 
ing regime only (white regions in Fig. |l|, top). Since an 
IR wave is a fundamental eigenstate of the self-induced 
waveguide, it can be demonstrated that the conditions 
of the Vakhitov-Kolokolov (VK) stability theorem |Tl| ] 
are satisfied, and the IR states are unstable if and only 
if dP/dui < 0. The critical point dP/duj = corre- 



sponds to a "collision" between two internal modes of 
the localized wave at the origin, as illustrated in Fig. |l| 
(middle). For the BR states, the VK criterion provides 
only a necessary condition for stability, since the higher- 
order localized states can also exhibit oscillatory insta- 
bilities. Indeed, we notice that in the linear limit there 
always exists an internal mode corresponding to a reso- 
nant coupling between the BR and IR band-gaps, since 
Y (^J^^ — lu^^^J — 0. We perform extensive numerical 

calculations and find that this mode leads to an oscilla- 
tory instability of BR waves when the value (lo — Re 7) 
moves outside the band gap; it happens when the inten- 
sity exceeds a threshold value (see Figs. |l| and ||). 




FIG. 4. Power vs. frequency, and (a,b) wave profiles for 
a — —5, (3 = —1. Notations are the same as in Fig. [I]. 

Next, we study the linearly repulsive defect (a < 
0) with self-focusing nonlinearity In this case, the 
waveguiding-type states bifurcate from the band-gap 
edges [with Go = ((uJb) = 0] and the stability criteria 
are the same as in the case a > discussed above, see 
Fig. |^. The major difference from the previous case is 
the existence of the anti-waveguiding states at small in- 
tensities, Jo < \ct/f3\. We find that such waves are stable 
in the lowest photonic band gap [such as one shown in 
Fig. [|(b)], while in higher band gaps they exhibit oscil- 
latory instabilities due to a resonance with the higher- 
frequency states. Different types of the localized waves 
and their stability are summarized in Fig. |§|. 

Self-defocusing nonlinearity. Let us now consider the 
case a < and (3 = — 1, when the total response of the 
defect is negative for any intensity, i.e. G(I) < 0. In 
this case, the localized waves can only exist in the anti- 
waveguiding regime (see Fig. ^j, top), the nonlinear waves 
continue the linear impurity states, and their frequency 
decreases at higher powers. We find that the localized 
states corresponding to the lowest BR gap, such as that 
shown in Fig. |](b), are always stable. However, in other 
band gaps the localized states [e.g., the one shown in 
Fig. |](a)] can exhibit oscillatory instability due to the 
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physical mechanism identified in other cases — a reso- 
nance with the lower band gaps. 

Finally, we study the case of a linearly attractive de- 
fect (a > 0) possessing a self-defocusing nonlinearity 
(f3 = —1). At small intensities, the defect is attrac- 
tive, and it can support localized waves in a waveguiding 
regime, in both IR and BR gaps, as shown in Fig. ||. 
Performing the stability analysis, we find that IR waves 
are always stable, while oscillatory instabilities appear for 
higher band-gap states. Since at higher intensities, i.e. 
for I > \a/fl\, the defect response changes its sign, a new 
type of localized waves can exist in the anti-waveguiding 
dispersion region, bifurcating from the band-gap edge. 
It is possible to demonstrate that for such waves the VK 
stability criterion becomes inverted, i.e. the waves are un- 
stable if dP/du> > 0, as shown in Fig. |S| (top). This hap- 
pens because the signs of both the nonlinear response and 
effective dispersion are altered compared to the IR waves 
supported by a self- focusing defect. All higher band-gap 
states can exhibit oscillatory instabilities. 
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FIG. 5. Power vs. frequency, and (a-c) wave profiles for 
a = 5, fH — — 1. Notations are the same as in Fig. |l| 

Concluding remarks. We have analyzed the exis- 
tence and stability of nonlinear localized waves in one- 
dimensional periodic structures considering the simplest 
nonlinear generalization of the Kronig-Penney model 
with a nonlinear impurity for which both the nonlinear 
wave classification and stability analysis can be carried 
out in a complete and systematic way. Taking into ac- 
count many common features of nonlinear guided waves 
and impurity modes in stratified and disordered media, 
on one side, and the self-trapped states and solitary waves 
in homogeneous nonlinear media, on the other side, we 
expect that many of the results described above will be 
found in more realistic physical models of nonlinear peri- 
odic media. Such cases include the electron self-trapping 
and locking states in cuprates and semiconductor super- 
lattic es |12| ; nonlinear guided waves in optical superlat- 
tices [13[; impurity modes in photonic band-gap mate- 
rials 14 1, magneto-optical periodic structures Jig], and 



photonic crystal fibers |]16); coherent matter waves in 
optical lattices jlTj, etc. In particular, the stability of 
nonlinear localized impurity modes is a crucial issue for 
creating tunable band-gap materials where gaps could be 
controlled by changing the input light intensity. 

We are indebted to O. Bang and C. M. Soukoulis for 
useful collaboration. The work has been partially sup- 
ported by the Planning and Performance Fund of the 
Institute of Advanced Studies. 
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